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Abstract. We present a theoretical model to describe the dynamics of Bose- 
Einstein condensates in anharmonic trapping potentials. To first approximation 
the center-of-mass motion is separated from the internal condensate dynamics 
and the problem is reduced to the well known scaling solutions for the Thomas- 
Fermi radii. We discuss the validity of this approach and analyze the model for 
an anharmonic waveguide geometry which was recently realized in an experiment 

m 

It is well known, that for a cloud of particles that is moving in a parabolic potential 
the center-of-mass motion is completely decoupled from the internal excitations [21 El • 
This holds in the presence of interaction between the particles and therefore also 
applies to Bose-Einstein condensates. Most of the experimental work with Bose- 
Einstein condensates has been done in parabolic potentials. The static properties of 
a condensate, see e. g. Ref. ^E], as well as the dynamical properties 013 have been 
investigated and were found to be in good agreement with the theoretical predictions 
|51I^ I1U[I11| . Higher order contributions to the trapping potential were either avoided 
or turned out to be negligible due to the large scale magnetic field generating elements. 
Anharmonic configurations were theoretically considered for a Gaussian 12 and a box- 
like potential and were experimentally realized in the case of a magnetic mirror 
|14| . Current progress in generating Bose-Einstein condensates in magnetic microtraps 
| 15l 115) in which the trapping potential in general differs from the trivial parabolic 
shape raises the questions how the dynamics of the trapped condensate are affected 
by the anharmonicity of the potential. In a recent experiment, the dynamics of a 
Bose-Einstein condensate in an anharmonic waveguide were studied in detail 1 and 
the theory in this paper has been developed to describe the experimental data. 

In anharmonic traps the center-of-mass motion of a condensate is coupled to 
the internal dynamics. In the rest frame of the condensate, the local potential 
curvature changes during the center-of-mass oscillation and collective modes of the 
condensate are excited. We present a theoretical model which allows us to describe 
the experimental data to high accuracy. The model is based on existing exact solutions 
of the Gross-Pitaevskii equation. Kagan et al. |E| and Castin and Dum developed 
an analytical solution of the time-dependent Gross-Pitaevskii equation for anisotropic 
time-dependent harmonic potentials in the Thomas-Fermi approximation. The time 
evolution of the Thomas-Fermi radii is thereby described by three ordinary nonlinear 
differential equations. To apply this solution we take advantage of the separation of 
the center-of-mass motion, following the lines of Japha and Band ^HI- After deducing 
the equations of motion we discuss the validity of this approach. We then apply the 
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theory to an anharmonic waveguide as it was studied in the experiment. We conclude 
by introducing the corresponding Hamiltonian formahsm. 

We start from the time-dependent Gross-Pitaevskii equation (GPE) 

t2 
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V2 + C/(x)+g|7^(x,t)|^ 



V^(x, t). 



(1) 



We introduce a vector R(t) for an arbitrary displacement of the condensate wave 
function and look for solutions of the form 



(2) 



(3{t) denotes an additional global phase factor and the gradient of the phase is given 
by the velocity R(t) of the displaced wave function 

TO, 



V5(x,i) = ^R. 



(3) 



Without loss of generality, the potential can be expressed in a Taylor series about R 



[/(x) = t/(R) + (x'V){/(R) + i(x'V)2c/(R) 



In equation Q we have used the abbreviation 



(4) 



(x'v)"= <A.■■■<AA.■■■^^^■ 

The shifted coordinates are given by 

Xj^ — 

and the spatial derivatives are identical in both coordinate systems 
d _ d 
dx[ dxi 

Inserting the wave function and the potential into the Gross-Pitaevskii equation 
CQ), we separate the result into three equations. The global phase I3{t) is given by 

ft 



(5) 

(6) 
(7) 



dt' 



-R(W)-^R(i')'--c/(R(0) 

2 TO 



For the displacement vector R we find 

toR = -V?7(R). 
Summing up all others terms leads to 



(8) 
(9) 



in^Jo{^',t) = - |-VVo(x',<)+5lV'o(x',i)|Vo(x',t) 



dt 



2m 



n=2 



(10) 



Note, that the sum in equation Hl()|l starts from n = 2 and the equation is written in 
the displaced coordinates x-. 

We now assume that higher order contributions in the sum of equation (|10|l 
are small compared to the quadratic term, which means that the potential can be 
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considered as locally quadratic. Then, taking only the leading term in the sum of 
equation (|10|) . one recognizes the GPE in the rest frame of the condensate 



d 



2m 



3 I 

1=1 



\t)x';'+g\M^',t)\' 



The time-dependent harmonic potential is given by 



1 



1,2,3, 



(11) 

(12) 



provided that the potential separates in the three coordinates and the mixed 
derivatives vanish. The solution of equation (|11|) within the Thomas-Fcrmi 
approximation is given in Refs. jl7l lll|. It consists of three scaling equations for 

the Thomas-Fermi radii of the condensate: 

,2 



A,; = 



-Lot{t)X,, z = l,2,3. 



(13) 



AiAiA2A3 

The Xi are scaling factors for the condensate radii ri(t) — A^rjo and are the Thomas- 
Fermi radii aX t — 0. The initial conditions are Ai(0) — 1 and Ai(0) = 0. The solution 
in the rest frame is given by |llj 




(14) 



X exp I — 



Equation H13I) is also valid for complex oscillation frequencies which corresponds 
to negative curvatures |19j . Thus, the dynamics of the condensate are completely 
determined by the set of ordinary differential equations ©, H12() and H13() . To get 
the complete condensate wave function in the laboratory frame one has to insert the 
shifted coordinates © into H14|) . Together with Q and (j2Jl this yields 



V'(x, t) 




X exp 



(15) 



The first line of p5|l is the square root of the density distribution of the condensate. 
As it is symmetric about x = R, the vector R coincides with the center of mass and 
therefore equation describes the motion of the center of mass. The second line of 
P5|l contains the spatially varying phase factors, arising from the internal dynamics 
(first two terms) and the center of mass motion (third term) . Their gradient defines 
the velocity field of the condensate. In the last line of H15() . all global phase factors 
are collected. 

We discuss briefly the consequences of the theoretical model. For a parabolic, 
linear or vanishing potential, the right hand side of expression H12|l is constant or 
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zero and the solution (|15|l is exact within the Thomas-Fermi approximation. The CM 
motion is completely decoupled from the internal degrees of freedom. For anharmonic 
potentials, equation l|12|) establishes the coupling between the external and internal 
dynamics. Equation Hll|l was derived under the assumption that terms of third 
and higher order can be neglected in the sum of equation H10() . This is a valid 
approximation as long as the potential curvature varies only slightly over the extension 
of the condensate. If the condition 

^f/(R)ro « ^C/(R) (16) 

is fulfilled the potential can be considered as locally quadratic. The inequality ((TB|l is 
violated when the quadratic part of the potential is zero or becomes negative. In this 
case the approximation remains valid if the contributions of the potential energy in 
pU|l to the total energy are much smaller than the chemical potential and the kinetic 
energy. 

We now apply the theory to an anharmonic waveguide. This geometry is the 
object of an actual study of the nonlinear condensate dynamics The potential 
of the waveguide is anharmonic in the z-direction and cylindrically symmetric in the 
yO"direction: 

= + az^ + + C2^, (17) 

where we have expanded the potential up to the fourth order. With the initial 
conditions -Rp(O) — 0, Rz{Q) — Ro and R(0) = the solution of (|5Jl can be written as 
a Fourier series 

oo 

Rz{t) — ^a„cos(nwoO (18) 

n=l 

Rp{t) = (19) 

Wo denotes the fundamental frequency, which depends on the amplitude of the center- 
of-mass oscillation. In the case of time of flight measurements (duration t), the 
amplitudes and phases of the Fourier components in the series ()18|l change: 



a'„=a„v/l+^ (20) 
4>[-^ = 4>n + arctanwnt. (21) 

Inserting ()18|l and H17() into the coupling equation H12() one obtains 
io^{tf = — [2a + 6bRz{t) + l2cRl{t)) 



2a + 6b a„ cos (nujot) + 12c ( ^ ^ cos (nujot) 

n=l 



1 



(2a + 6cai + 66ai cos{uJot) + %ca\ cos(2ajo0) • (22) 



In the second step we have made the approximation, that ai ^ a„ for n > 1, which 
means, that the dominant Fourier component is at the fundamental frequency ojq. In 
the rest frame, the condensate is driven in its axial direction with ojq. The driving 
amplitude is directly proportional to the amplitude of the center-of-mass motion ai. 
Due to the contribution of the fourth order in the potential Ijl?! . the condensate is 
additionally driven at twice the fundamental frequency with a quadratic dependence 
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on ai- Beside these two frequencies, the (m=0) low-lying quadrupole resonance is 
offresonantly excited For small amplitudes, the oscillation frequency in the rest 
frame becomes constant and equal to the frequency of the center-of-mass motion 
{iOz = ujQ — yjlajni). For strong excitation, the nonlinearity of the equations of 
motion for the internal condensates dynamics H13() produces a nonlinear coupling of 
the collective excitations. A comparison with actual experimental data shows an 
excellent agreement between theory and experiment 

Kagan et al. JJj pointed out, that the equations of motion for the internal 
condensate dynamics H13f) can also be derived from a classical Hamiltonian. For the 
free motion in a trapping potential the total energy of the condensate is conserved 
and can be calculated for the wave function (|15|l 

H = (Skin) + (£;pot) + (Sint) 

+ im^i?2 + C/(R),. (23) 

i 

The first three terms in H23|l are the contributions from the internal dynamics and the 
last two terms that of the center-of-mass motion. Taking equation H23|l as a classical 
Hamiltonian with the canonical variables 

= {i?i,i?2,i?3,Ai,A2,A3} (24) 
and the conjugated momenta 

Pi = {mRi,mR2,mR:i, -;^^>^2, t^^As}. (25) 

'^10 '^20 ''^30 

one can evaluate the Hamilton equations of motion 
dH 
opi 
dH 

We obtain the equations H13|l and the definition of the time-varying oscillation 
frequencies (|12|l . For the centcr-of-mass motion we find a slightly modified equation 
of motion 

-8 1 (9^ 

The second term in (|26|) is a small correction to Q which accounts for the influences of 
the internal dynamics to the center-of-mass motion. The correction is proportional to 
the square of the condensate radii, whereas the internal dynamics H13|) are independent 
of the extension of the condensate. The modified equation of motion (|26(l is useful to 
test the numerical stability of the calculation although its influence on the amplitude 
of the center-of-mass oscillation is negligible. 

In conclusion, we have presented an approximate solution to the time-dependent 
Gross-Pitaevskii equation for anharmonic potentials. We have shown, that the internal 
dynamics of the condensate are coupled to the external center-of-mass motion. For 
an anharmonic waveguide, we have analyzed the internal dynamics in detail and good 
agreement with the experimental data is found P . We have also introduced a classical 
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Hamiltonian from which the equations of motions can be derived and have discussed 
the validity of the approximations. 
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